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ABSTRACT 

In this paper we present the results obtained with a monitoring programme (23 days 
long) performed with Swift-XKT on the local Seyfert galaxy Mrk 915. The light-curve 
analysis shows a significant count rate variation (about a factor of 2-3) on a time- 
scale of a few days, while the X-ray colours show a change in the spectral curvature 
below 2 keV and the presence of two main spectral states. From the spectral analysis 
we find that the observed variations can be explained by the change of the intrinsic 
nuclear power (about a factor of 1.5) coupled with a change of the properties of an 
ionized absorber. The quality of the data prevents us from firmly establishing if the 
spectral variation is due to a change in the ionization state and/or in the covering 
factor of the absorbing medium. The latter scenario would imply a clumpy structure 
of the ionized medium. By combining the information provided by the light curve and 
the spectral analyses, we can derive some constraints on the location of the absorber 
under the hypotheses of either homogeneous or clumpy medium. In both cases, we find 
that the absorber should be located inside the outer edge of an extended torus and, 
in particular, under the clumpy hypothesis, it should be located near, or just outside, 
to the broad emission line region. 
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1 INTRODUCTION 

There is now a general consensus that active galactic nuclei 
(AGN) are powered by accretion of matter on to a super- 
massive (>10® Mq) bl ack hole (SMBH) . According to the 
unified model of AGN J Antonuccil 1 1 9931 1 . an obscuring op¬ 
tically thick medium composed by dust and gas arranged 
in a torus-like geometry is present around the nuclear en¬ 
gine. However, the structure, size and composition of this 
circumnuclear medium are still matter of debate and are 
topics of several studies carried out at different wavelengths 
iBianchi, Maiolino. fc Risaliti|[20l3 l. Important constraints 
on the physical properties of the circumnuclear medium have 
been recently provided by the study of the absorption vari- 
ability, which is alm ost ubiquitous in bright absorbed AGN 
llRisaliti et al.ll^02ll . X-ray absorbing column density (Nh) 
has been observed to vary by a factor of 10 or more over a 
few years. The observed Ah variability on this scale ruled 
out the simplest physical configuration of a homogeneous 
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absorber, giving upper limits on the distance of the latter 
by the central SMBH. Subsequent X-ray observational cam¬ 
paigns detected significant variability on a handful of AGN 
on very short time-scale (day or week) giving the possibil¬ 
ity to investigate the X-ray properties of the circumnuclear 
medium down to sub-parsec scale. The emerging picture is 
that multiple neutral and ionized absorbers co-exist around 
the central SMBH, located at different distances from it. So 
far the best characterization of the physical parameters of 
the X-ray absorbers has been possible for the AGN in NGC 
1365, which has been monitored several times in the last 
few years with Chandra, XMM-Newton and Suzaku. This 
source changes from Compton-thick (AhIsIO^'^ cm“^) to 
Compt on-thin (Ah^IO^^ cm~^) state on ti me-scales from 
week s ( Risaliti et ahlf^OSll to hours or days llRisaliti et al.l 
l2007l . l2009h . One of the main results is that the X-ray ab¬ 
sorber has the physical properties typical of the clouds re¬ 
sponsible for the emission of the broad lines in the opti- 
cal/UV spectrum, i.e. the broad-line emission region (BLR), 
located at a distance of hundreds of gravitational radii from 
the central SMBH. Recently it has been proposed that the 
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long-term variability of the X-ray absorber in NGC 1365 
could be due to the variation of the ionized X-ray wind, 
that is responding to the changes in the accretion rate 
iConnollv. McHardv. fc Dwellvll20Ml . High spectral resolu¬ 
tion data show that this wind is located within the variable 
soft X-ray absorber and is composed by two (or even three) 
zones with different ionization levels. In particular, the low¬ 
est ionization zone of this wind could be responsible for the 
absorp tion variability tha t is observed in this source below 
2 keV llBraito et al.H2oi3i . 

By exploiting th e Swift’K-ray Telescope (hereafter XRT; 
[Burrows et al.ll2005ll archive, we recently started a project 
aimed at studying, on a larger statistical basis, the physi¬ 
cal properties of the X-ray absorbers in AGN. We started 
from the 70 months all- sky survey Swift-HAT catalogue 
([Baumgartner et al.ll201^ in the 20-150 keV and the 66- 
month Palermo BAT Gatalogu^il (see also [Cusumano et al] 
[2010lj . We considered only those sources with a secure op¬ 
tical AGN counterpart and with XRT observations already 
available before the end of 2012. For each source, we adopted 
a daily binning of the XRT data and we searched for statisti¬ 
cally significant variations of the X-ray colours on differen t 
time-scales (from months to years, see [Ballo et al.[ [2015li . 
Following this method, we selected five new candidates for 
variable absorbers. One of the most interesting is the local 
Seyfert galaxy Mrk 915. The archival XRT data used to se¬ 
lect this source unveiled a dramatic count rate and ‘X-ray 
colour’ variability on a minimum time-scale of 1 month. The 
0.3-10 keV count rate increases by a factor of 3 and, at the 
same time, the spectral shape becomes dramatically softer, 
showing a variability of a factor of 2.4 and 3 in the [2-10 
keV]/[0.3-2 keV] and [4-10 keV] /[0.3-4 keV ] count rate ra¬ 
tios, respectively. As discussed in iBallo et al.l l[2015ll . this be¬ 
haviour suggests that a change in the properties of the X-ray 
obscuring circumnuclear medium (i.e. column density, cov¬ 
ering factor or ionization state) probably occurred between 
the archival XRT observations. However, the low statistical 
quality of these data prevented us from performing a proper 
spectral analysis and to discern among the possible causes 
of the observed variation. With the goal of obtaining higher 
quality data on a shorter time-scale, we performed a daily 
XRT monitoring on this source during a period of ~3 weeks. 
This paper aims at providing a first characterization of the 
putative variable absorber. 

The paper is structured in the following way: the main 
optical properties of Mrk 915 are reported in Section 2, while 
the analysis of the XRT data obtained during the daily mon¬ 
itoring are presented in Section 3 (light-curve analysis and 
spectral analysis. Sections 3.1 and 3.2, respectively). The 
main results obtained are discussed in Section 4 and the 
conclusions are summarized in Section 5. 

Throughout the paper we assume a flat AGDM cosmol¬ 
ogy with Ho=71 km s“^ Mpc“^, Ha=0.7 and Hm=0.3. 

2 MRK 915 

Mrk 915 is a local (j;=0.024) spiral galaxy (PA~166°) 
with evident dust lane structures seen crossing the cen- 

^ http://bat.ifc.inaf.it/bat_catalog_web/66m_bat_catalog.html 


Table 1. Mrk 915: XRT monitoring observation log. The obser¬ 
vations are ordered on the basis of the observation start date. 


Obs. ID 

Obs. start date 

Net counts 
[0.3-10 keV] 

Net exp. time 

[s] 

00035169003 

2014-09-10 

638 

7367 

00035169004 

2014-09-11 

621 

7060 

00035169005 

2014-09-15 

344 

4453 

00035169006 

2014-09-17 

1359 

13130 

00035169007 

2014-09-18 

437 

3629 

00035169008 

2014-09-21 

2563 

13190 

00035169009 

2014-09-23 

693 

4105 

00035169013 

2014-09-27 

160 

1171 

00035169014 

2014-09-27 

61 

719 

00035169011 

2014-09-29 

935 

10280 

00035169015 

2014-10-02 

993 

5369 

00035169016 

2014-10-02 

1216 

5314 


tral s ource in its WFPC2 image l|Malkan. Goriian. fc Tarnl 
[l998ll and hosting a n SMBH of Mbh=(0.6-1.8) x 10® Mq 
i Bennert et al.[[2006tl. The v alue of the axis ratio (0.54 and 
0.83 reported bv [Keell[ 19961 and [Munoz Marfn et al.[[2007[ . 
respectively) indicates that the disc of the galaxy is inclined 
with respect to the line of sight by about 35°-57°. Assum¬ 
ing both galaxy’s disc and AGN’s torus to be co-planar and 
an half-opening angle of the torus of about 60°, our line of 
sight intercepts or grazes the low-density outer regions of 
the torus. 

From the optical point of view, Mrk 915 has been spec¬ 
troscopically classified as an intermediate SeyferlQ galaxy 
showing significant spectral variations. In particular, the 
broad components of the Ha and H/I emission lines are 
clearly detected in some observations, while they are com¬ 
pletely absent in others. As a consequence, the source has 
been classified in different ways, from Seyfert 1.5 to Seyfert 
1.9, depending on the presence/intensity of the broad com¬ 
ponents at the time of the observatio ns ([Goodrichl [ 19951 : 
[Bennert et al.[ [20061 : [Trippe et al.[ [ 2 OI 0 II. The origitr of this 
spectral variability is still unclear. [Goodrichl ([l995^ tenta¬ 
tively suggested that the observed variation was due to a 
change in reddening (AEfl_v^0.53 mag, derived by com¬ 
paring the Ha fluxes in the two observations) produced by 
dusty clouds passing close, but outside, the bulk of the BLR. 
Alternatively, the observed spectral variability could be due 
either to a change in the nuclear photoionizing continuum 
or to a combination of continuum and reddening variation. 
At the moment, it is not possible to d iscriminate among the 
various scenarios ([Trippe et al.[[201^ . 

Independently of the origin of this spectral variability, 
a reddening of the order of Eb-v~0.3-0.5 mag seems to be 
always present and it affects both the kpc-distance regions, 
i.e. the narrow-line emission regions (NLR), and the inner¬ 
most nucleus. This reddening is most probably attributed to 
the presence of dust l anes as seen crossing the ce ntral source 
in its WFPC2 image ([Munoz Marfn et al.[[200^ 1. 


^ A Seyfert galaxy with broad emission line com ponents weake r 
than those usually observed in Seyfert 1 galaxies liWhittlal 19921) . 
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Figure 1. XRT light curve in the total 0.3-10 keV band, obtained by binning the data per snapshot (small dots, grey in the electronic 
version) and per observation (big dots, red in the electronic version). With respect to the observed net counts reported in Table[T] data 
points are corrected for technical issues (i.e. bad p ixels/column s , field of view effects, pile-up and source counts landing outside the 
extraction region) following the receipts discussed bv lEvans et al.l l2^^7^ and lEvans et all ll20ni3ll . Numbers mark the one/two last digital 
numbers of the relevant observation ID (see Table (TJ. Error bars mark Icr uncertainties. The light curve is divided in time intervals 
marked by different boxes and capital letters (see Section 3.2). We note that the only point in the B time range with a count rate lower 
than 0.2 has a very short exposure time (~180 s). Although it is above the minimum exposure time considered here (i.e. 150 s), this 
point has a signal—to-noise ratio lower than 5. 




Figure 2. Left panels, from top to bottom: 0.3-2 keV and 2-10 keV count rates and [2-10 keV]/[0.3-2 keV] ratio. Right panels, from 
top to bottom; 0.3-4 keV and 4-10 keV count rates and [4-10 keV]/[0.3-4 keV] ratio. Dashed lines represent the ratios normalized to 
the first observation. Symbols are the same as those used in Fig. [T] and data points are corrected for the same technical issues. 


3 DAILY/WEEKLY X-RAY MONITORING 

The data presented here are relative to a monitoring pro¬ 
gramme composed by 12 new XRT observations awarded 
dnring fall 2014 {Swift Cycle-10 and ToO observations, P.I. 
Severgnini). These observations, covering ~23 days, for a 
total on source exposnre time of ~76 ks, were performed 
with XRT in the standard PC-mode. The observation log 
is reported in Table [T] The sonrce appears point-like in the 


XRT image and we do not find any significant evidence of 
pile-up. 

3.1 Light curves 

Each light curve, described below, has been created by us¬ 
ing the Swift-XRT data products generatoi0 developed by 

® http://www.swift.ac.uk/user_objects/ 
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Energy (keV) 

Figure 3. XRT spectra relevant to the time binning defined in 
Fig.m In the electronic version, the spectra are colour-coded as 
the boxes in Fig. [T] 

the UK Swift Science Data Centre and based on the Swift 
software and ftools guide. The effects of the damage to the 
CCD and automatic readout-mode switching are appropri- 
at ely handled; details on t he light-curve code can be found 
in lEvans et ahl (|2007l ~) and lEvans et al.l (l200fll 'l. 

Figs [T] and [2] show the XRT light curves in five different 
energy bands (0.3-10, 0.3-2, 0.3-4, 2-10 and 4-10 keV). We 
choose to bin the data per spacecraft orbit (i.e snapshot, 
grey small dots in Figs [T] and [2)). We considered only those 
snapshots longer than 150 seconds, which guarantee a signal- 
to-noise ratio on each point equal or greater than 5. In each 
panel, we plot also the light curves obtained by binning the 
data per observation (i.e. obs. ID, big red dots in Figs[T]and 
El. Over the length of the observing monitoring programme 
(~23 days), the light curves are not constant at 99.9 per 
cent confidence level (x^ test). 

During the first three observations, spanning ~3.6xl0® 
s (i.e. 4.2 days), the snapshot 0.3-10 keV count rates are ran¬ 
domly distributed around ~0.12 counts s“^. Later, a clear 
increase in the count rate is visible, reaching its maximum 
in about 3.5x10® s (i.e. ^4 days). In these four days, the 
count rate increased by a factor ~2-3 when we consider the 
obs. ID time binning. About 2 days later, the count rate de¬ 
creased by the same amount and then immediately increased 
again, reaching the highest intensities observed during this 
monitoring (~0.3 counts s“^). 

The increasing/decreasing factor is not the same in all 
bands: in the softer bands (0.3-2 keV and 0.3-4 keV) the 
count rate variation is higher with respect to that in the 
complementary harder bands (2-10 keV and 4-10 keV, re¬ 
spectively). This is visible also inspecting the flux ratios ([2- 
10 keV/0.3-2 keV] and [4-10 keV/0.3-4 keV]) reported in 
the bottom panels of Fig. El In these panels, the dashed lines 
are normalized to the first observation. This behaviour sug¬ 
gests that, besides the possible flux variation of the primary 
emission, there could be another mechanism producing dif¬ 
ferent intensity variations in different energy bands. This 
mechanism should affect mainly the softer energy ranges. 

3.2 X-ray spectral analysis 

The aim of this section is to investigate, from a spectral point 
of view, the changes observed in the light curves of Mrk 915 
during the X-ray monitoring. To this end, following the to¬ 


tal (0.3-10 keV) light curve, we divided the XRT data in 
four main time bins on the basis of the average count rate 
value of each observation. They are shown in Fig. [H where 
we mark each bin with boxes and capital letters. The A and 
C bins are composed by observations with a net count rate 
lower than -^0.2 counts s“^, while B and D correspond to 
observations with a net count rate higher than ~0.2 counts 
s“^. For each bin, we extracted the spectrum by using cir¬ 
cular regions centred on the X-ray source position with a 
radius -^50 arcsec (i.e. ~20 pixel, which corr esponds to an 
Encir cled Energy Fraction of ~90 per cent; iMoretti et al.l 
l2005h . The background spectra have been extracted from 
source-free circular regions close to the object with an area 
about 9 times larger. Spectral reduction was performed us¬ 
ing the standard software (headas software, v6.15 and the 
most updated CALDB versiorQ) and following the proce¬ 
dures described in the instrument user guid^B- The spectra 
have been binned in order to have at least 20 counts per 
energy channel (see Fig. El and analysed using the xspec 
12.8.1 package. 

As evident from Fig. El tke data have two by two (A 
with C and B with D) very similar spectral shape, confirming 
the presence of two different spectral states of the source, as 
already highlighted by the light curve and X-ray colours. A 
visual inspection of the figure is sufficient to notice that most 
of the spectral variation among these two states is below 2- 
3 keV. In order to improve the statistics, we combined the 
data sets with similar spectral shape (A with C and B with 
D) by obtaining a low and a high-state spectrum (hereafter 
called AC and BD, respectively). 

We also considered the Swift-BAT spectrum of Mrk 
915 retrieved by t he 70 month Swift-BAT Catalogu43 
llBaumgartner et al.ll20R^ . This spectrum provides the av¬ 
erage spectral properties of the source above ~15 keV. Due 
to the significant variability of Mrk 915, we did not fit the 
XRT and BAT spectra simultaneously. We used the BAT 
spectrum to derive the average value of the spectral slope of 
the primary AGN emission. By htting the BAT data with a 
simple power-law component, we found a best-fitting value 
of r=1.7±0.2. In the spectral analysis of the XRT data (de¬ 
scribed in details below), we hrst htted simultaneously the 
two main states keeping F tied. Also in this case we found: 
F = 1.7 ±0.2. The same value and uncertainties are obtained 
when F is allowed to vary between the states. Thus, consid¬ 
ering the statistic of the present data, we decided to keep it 
fixed to the best-fitting value while investigating the origin 
of the variability. 

As a starting point, we assumed the presence of three 
absorbing components. The first one {phabs component in 
xspec) refers to th e Galactic hydrogen co lumn density along 
the line of sight dKalberla et al.l l2005h . The second one 
{zphabs component in xspec) accounts for the absorption 
due to the dust lane structures at the redshift of the source. 
Indeed, as discussed in Section 2, the primary emission 
from the Mrk 915 nucleus intercepts dust lane structures 
that produce an opt i cal extinction of E b-v= 0.3-0.5 mag 
jBennert et al.1120061 : iTrippe et al.1 12 OI 0 I ') . The presence of 


http://heasarc.gsfc.nasa.gov/FTP/caldb 
® http://heasarc.nasa.gov/docs/swift/analysis/documentation 
® http://swift.gsfc.nasa.gov/results/bs70mon/ 
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Table 2. Best-fitting values obtained by applying to the low- and high-state spectra of Mrk 915 the 
following model: phabsx{zphabsxzxipcf xzpowerlw) with F fixed to 1.7 (x^/d.o.f.=362.99/403). 


State 

neutral 

ionized 

CF= 

log(^)'^ 

Ae 

F(2-10keV)'^ 

7^{2-10keV)® 

AC 

0.08t°°3 

1.86±0.20 

0.96±0.03 

0.68±0.13 

1.98±0.01 

0.67 

0.98 

BD 

0.08'* 

2.09±0.35 


-*-•-*-^-0.24 

3.06±0.01 

1.07 

1.52 


Notes: The uncertainties, reported at the 90 per cent confidence level for one parameter of interest, 
have been estimated by freezing the A^Hneutral to the best-fitting value. 

“ Column density of the neutral absorber associated with dust lane structures, in units of 10^^ cm~^. 
This parameter is forced to be equal in both states. ^ Column density of the ionized absorber, in 
units of 10^^ cm“^. ^ Covering factor of the neutral absorber. ^ Logarithm value of the ionization 
parameter (see Section 3.2). ® Normalization of the power-law component, in units of 10“^ Pho¬ 
tons s~^ cm~^ keV“^. ^ Flux in the 2-10 keV energy range, corrected for the Galactic absorption, 
in units of 10“^^ erg cm“^ s~^. ^ Intrinsic 2—10 keV luminosity, in units of 10"^^ erg s“^. ^ This 
parameter was kept tied in the two states. 





I _I. . .I_I_I_ 

' 0.5 1 2 5 

Energy (keV) 


Figure 4. Residuals, plotted in terms of sigmas with er¬ 
ror bars of size one, obtained by fitting the low- (AC) and 
high-state (BD) spectra of Mrk 915 with the following model: 
phabsx {zphabsx zpcfabsx zpowerlw). 


these structures is expected to affect also the X-ray emission 
with a column density of the order of 10^^ cm“^, assuming 
a Ga lactic dust-to-gas ratio of Eg- i//iVH=1.7xlO"2^ mag 
cm^ iBohlin, Savage, fc Drake! 1 19781 ). In the fitting proce¬ 
dure, we allowed this component to vary in the range of 
0.5-5x10^^ cm“^. The third component is a neutral ab¬ 
sorber partially covering the central source {zpcfabs model 
in XSPEC), associated with material on the torus distance 
scale. The full parametrization adopted is the following one: 
phabsx (zphabsx zpcfabsx zpowerlw). We allowed to vary the 
intrinsic power-law normalization and both the intrinsic col¬ 
umn density and the covering factor of the partial cover¬ 
ing absorber. The residuals between the data and the best¬ 
fitting model for the two states are shown in Fig. [4] 

Although this model can account for most of the data, 
it leaves evident and systematical residuals in the softer 
part of the spectra: by considering the data in the 0.3-1 
keV range we found x^/ d.o.f.=55.42/37 (corresponding to 
a null hypothesis probability of -^3 per cent). Thus, this 
model can not be considered a good representation of the 
global spectral properties of the source. In particular, the 
shape of the residuals in Fig. |4] suggests the presence of 
an absorption trough at ~0.78 keV (observed frame); at 
the redshift of the source, this feature could be associated 
with Fe XVII-XVIII, typically observed in nearby bright AGN 


t' 0.01 

> 

CD 

CO 

E 10-3 

o 

CO 

c 

o 

o 

f 10-“ 

> 

^ 4 

2 
0 

-2 


Figure 5. Upper panel: The model including an ionized absorber 
partially covering the source (phabsx (zphabsx zxipcfx zpowerlw)) 
is overplotted on the low- (AC) and high-state (BD) unfolded 
spectra of Mrk 915. Lower panel: Relevant residuals, plotted in 
terms of sigmas with error bars of size one. 

with warm absorbe rs l|Crenshaw. Kraemer. fc Georgd [20031 : 
IPorauet et al.l[2004h . Prompted by these considerations, we 
replaced the neutral partial covering compo nent with an 
ioniz ed one, modelled in XSPEC with zxipcf l|R,eeves et al.l 
I 2 OO 8 I) . In this case, the full parametrization is the follow¬ 
ing: phabsx (zphabsx zxipcfx zpowerlw). As a first step, in 
the fitting procedure, we left the intrinsic Vh, the ioniza¬ 
tion parameter ^ and the covering factor (CF) of the ion¬ 
ized medium free to vary, along with the intrinsic power- 
law normalizations (Ap). The best-fitting values obtained in 
the presence of an ionized absorber are reported in Tabled 
where the uncertainties are given at the 90 per cent confi¬ 
dence level for one parameter of interest llAvni|[l976l) . The 
unfolded spectra and the ratio between data and best-fitting 


^ The ion ization parameter is defined a s ^[erg cm s“^] = 
Lion/nl?^ llTarter. Tucker, fc Salpetej Il969l ). where Lion is the 
ionization luminosity obtained by integrating the X-ray intrinsic 
luminosity between 0.013 and 13 keV, n is the average absorber 
number density [part/cm“®] of the illuminated slab, and R is the 
distance of the absorber from the central source. 
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Figure 6. 68 per cent, 90 per cent and 99 per cent confidence 
contours for the AC state and BC states for different parameters 
of interest obtained by using the model reported in Table |2] Top 
panel: ionization state versus the power-law normalization (Ap). 
Middle panel: Covering factor versus the power-law normalization 
(Ap). Bottom panel: Covering factor versus the ionization state. 

models are shown in Fig. O it is evident that this second 
model provides a better characterization of the spectra of 
Mrk 915, also in the softer energy range (x^/d.o.f.=22.42/35 
in the 0.3-1 keV range, corresponding to a null hypothesis 
probability of ~95 per cent). Also the residuals at E ~0.78 
keV are now well reproduced. We note that, with the present 
data, we cannot investigate the presence of a possible out¬ 
flow velocity. Thus, the redshift of the ionized absorber was 
fixed to the same 2 : of the source. 


The confidence contour plots of the joint errors for dif¬ 
ferent parameters of interest (log(/), CF and Ap), obtained 
with the same choices of fixed parameters as quoted in Ta¬ 
ble [ 2 I are shown in Fig. (6] The emerging scenario is that, 
in both states, the nucleus is obscured by a variable ion¬ 
ized medium with an A^h of the order of 2x10^^ cm“^ and 
a significant increase in the power-law normalization by a 
factor of -^1.5 from the low- (AC) to the high-state (BD) is 
also required (see Table[2|. We note that, this increase is re¬ 
quired also in the case of neutral partial covering absorber. 
In combination with an intrinsic fading/increasing of the 
power-law normalization, Fig. [6] suggests also a marginally 
evidence for a change of the ionization state and of the CF, 
albeit with low significance. In particular, when the source 
gets brighter the ionization state tentatively increases (90 
per cent confidence level, see Fig. [O] top panel) and the CF 
slightly decreases (68 per cent confidence level, see Fig. [6l 
middle panel). 

In the bottom panel of Fig. [HI we report the confidence 
contour plot of the joint errors of the CF versus the ioniza¬ 
tion parameter. This plot indicates that our spectra could be 
equally reproduced by assuming that only one of these two 
parameters is changing. Thus, we re-fitted the data twice, 
forcing the two spectral states to have: (1) the same ion¬ 
ization parameter, allowing the CF to be free to vary, and, 
(2) the same CF, allowing the ionization parameter to be 
free to vary. As expected, in both the cases we were able 
to reproduce our spectra with a statistical quality similar 
to that reported in Table [2] and showed in Fig. [5] In both 
cases, the intrinsic column density for the ionized medium 
(Ah~2x 10^^ cm“^) and the change in the power-law nor¬ 
malization (by a factor of 1.5) are consistent with the values 
previously obtained. 

In the case in which the two ionization parameters were 
tied together (log(5)=0.63±0.15), the spectral shape varia¬ 
tion could be completely ascribed to a change in the CF from 
-^0.97 in the low-spectral state to ~0.85 in the high-spectral 
state (x^/d.o.f.=365.33/404). In the case with the two cover¬ 
ing factors tied together (CF=0.95±0.03), the spectral shape 
variation could be completely ascribed to a change in the 
log(^) from ~0.77 in the low-spectral state to -^1.35 in the 
high-spectral state (x^/d.o.f.=363.5/404). 

Finally, we tested the picture discussed above fitting our 
final model, with both hypotheses of constant / or constant 
CF, to the single spectra showed in Fig.O(i.e. A, B, C, and 
D spectra spaced by about 5x10® s, see Fig. [T]). We found 
results consistent with those previously described. In spite 
of the lower statistics of the data with respect to the com¬ 
bined spectra, the variations in the ^ and CF parameters 
still remains statistically significant (68-90 per cent confi¬ 
dence level) as the change in the power-law normalizations 
(more than 99 per cent confidence level). 


4 DISCUSSION 

The spectral analysis performed on the XRT monitoring 
data of Mrk 915 unveiled the presence of an ionized absorber 
with Nh of the order of 2x10^^ cm“^ and allowed us to de¬ 
tect a variability of the intrinsic power of the central source 
of a factor of -^1.5 on a time-scale of a few days. The pos¬ 
sibility that the source could intrinsically decrease/increase 
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its power was already suggested from observations in the 
optical domain (see Section 2 for relevant references). 

As for the spectral variation observed in the X-ray do¬ 
main, we found that it could be due to a change in the 
ionization state of the circumnuclear medium and/or to a 
change of the covering of absorbing medium along the line 
of sight. Higher quality data are needed to understand if 
only one or both these parameters actually vary. If the vari¬ 
ation of the covering factor will be confirmed, this implies 
a clumpy structure of the absorber. Since our X-ray data 
can not firmly confirm or discard this hypothesis, in the fol¬ 
lowing we discuss some constraints on the location of the 
absorber by assuming both a homogeneous and a clumpy 
medium. 

Homogeneous medium. Under the hypothesis of homo¬ 
geneous medium, we can place an upper limit on the location 
of this absorber by using the relation between the intrin¬ 
sic continuum luminosity, the ionization and the density of 
the medium: ^ = Lion/nTJ^. Assuming that the thickness 
of the absorber is less than the distance from the central 
SMBH (AR/R< 1), this relation provides an upper limit on 
its distance: RKLion/NnC From the X-ray spectral analysis 
we derived: Ah~2x 10^^ cm“^ and Lion=(2.4-3.7) x 10"^®erg 
s“^ (low-high state, respectively; Lion is the intrinsic X- 
ray luminosity, corrected for absorption, and integrated be¬ 
tween 0.013 and 13 keV, see also footnote 7). By considering 
the minimum value of ^ obtained under the hypothesis of 
no covering factor variation, we obtain i?< (2-3) x 10^° cm, 
i.e. lower than (0.8-3)xl0^ i?g {Rg = GMBH/c^ is the grav¬ 
itational radius, where Mbh is the SMBH mass equal to 
0.6-1.8x10® Mm: iBennert et ^120061 '). This estimate corre¬ 
sponds to a location inside the outer edge of the torus. 

Clumpy medium. Under the hypothesis of a clumpy 
medium and assuming that the absorber is composed by 
spherical gas clouds m oving around the central source with 
Keplerian velocity (e.g. lRisa.liti et al.ll2002l . lElvis et al.l[2003 . 
iRisaliti et al.l 12003 1. the variability time-scale between the 
two maxima (B and D) detected in Fig. [T] and Fig. [21 
<t>=10® s, at the first order, is equal to the crossing time 
for a single cloud through the source. 

The distance Rc of the absorber from the central source 
is given by: Rc = GMbh / (Ds+Dcf where L is the 
variability time and He and Hs are the cloud and cen¬ 
tral source linear size, respectively. Under the hypothesis 
that the linear size of a single cloud is of the same order 
of that of the emitting central source, and taking into ac¬ 
count that Dc=Nii,c/Uc, we can re-write the distance Rc in 
the following way: Rc= 0.25 GMbb Nh,c~^, i.e., for 

A^h, 0=2x10^2 cm-^: 

Rc ~ 5 — 15 nf cm (1) 

The range of values reported in equation (1) is due to the 
uncertainty in the Mbh (0.6-1.8x10® Mq, see Section 2). 

A conceptually different estimate of Rc can be derived 
from the relation quoted above between the intrinsic contin¬ 
uum luminosity, the ionization parameter and the density 
of the absorber (^ = Lion/ucRc)- By considering an aver¬ 
age value of the ^ parameter of 10 erg cm s“^ and of the 
intrinsic ionization luminosity of Lion~3xl0^® erg s“^, we 
can derive: 

Rc = (3 X Uc^Y^^cm. (2) 


By combining the results reported in (1) and (2), we find 
that, under the hypothesis of absorbing clouds with lin¬ 
ear dimension similar to that of the central source, the 
clumpy medium should be located at a distance Rc = (1.4- 
1.7)xlO^’^ cm (i.e. ~10^Rg), with an average cloud density 
of nc~(l-1.7) X10® cm~® and an average linear dimension of 
Hc~Hs~(1-1.2)x 10^'‘ cm (i.e. '-7Rg). 

We repeated the estimates discussed above by taking 
into account that the linear size of the single cloud could in¬ 
deed be smaller (e.g. Hc~0.5 Hs) or greater (Hc>>Hs) with 
respect to that of the central source. We found that Rc could 
vary between ~10^^ cm and 2.2x10^^ cm with a density ric 
between 2.2x10® cm“® and 6x10^ cm~®, respectively. 

By considering the hlOOA luminosity estimated by 
iBennert et al.l (l2006h for this sourc e (ALa 1.68 x lO"^'^ 
erg s“^) and the relation found by iBentS (l2009l l between 
the BLR size and the 5100A luminosity, the BLR in Mrk 
915 is expected to be located at a distance from the cen¬ 
tre SMBH of the order of 10^®- 10^^ cm. If confirmed, the 
short time-scale variation of the CF would imply that the 
ionized absorber clouds should be located near, or just out¬ 
side, to the bulk of the BLR with a density between 6x10^ 
and 2.2x10® cm“®. This range of densities would indeed be 
consistent with the expected lower limit of the electron den - 
sity for the BLR (nc(BLR)>10® cm“®; lOsterbrockI ll98^ L 
We note that this s cenario is very similar to that proposed 
by iGoodrici] (Il995f l. see Section 2, to explain the optical 
variability of this source. 


5 SUMMARY AND CONCLUSIONS 

From the optical point of view, Mrk 915 is a spiral galaxy 
hosting an intermediate (type 1.5-1.9) AGN with evident 
dust lane structures crossing the central source. Significant 
variability (about a factor of 2.5-3) in the intensity of both 
the broad component of the Ha emission line and the un¬ 
derlying continuum has been observed and reported in the 
literature on time-scale of few years. The analyses carried 
out by different authors suggest that this change could be 
due to an increasing of the reddening along the line of sight, 
to a fading of the central source, or to a combination of 
them. Thanks to the unique capabilities of XRT on-board 
the Swift satellite, it was possible to investigate the daily 
time-scale variation of the X-ray absorber. 

As discussed in Section 1, Mrk 915 was selected as a 
good candidate for a variable absorber on the basis of the 
X-ray colour analysis on archival XRT data. The analysis 
of these observations (iBallo et al.ll2015ll . which caught the 
source in two different states, suggested that this could be 
associated with a variable absorber. 

With the aim at obtaining higher quality data on 
shorter time-scale, we performed a XRT monitoring cover¬ 
ing -^3 weeks and sampling days time-scale. During these 
new observations the source was caught again in two main 
different count rate and spectral states. We investigated the 
nature of the variable absorber testing two partial cover¬ 
ing models: in the first one the absorber is neutral, while 
in the second one it is ionized. In both cases, we added a 
neutral absorber (Ah~ 10^^ cm“^) associated with the dust 
lanes structures seen crossing the central source. We find 
that, while the first model (neutral absorber) is not able 


© XXXX RAS, MNRAS 000, HHS] 
























8 P. Severgnini et al. 


to reproduce the data below 1 keV, the second one (ion¬ 
ized absorber) provides a good characterization of the data 
in the full energy range considered here. Furthermore, we 
detect a significant variation of the intrinsic power of the 
nuclear source. The quality of the present data prevents us 
from firmly establishing if the observed spectral variations 
were due to a change in the ionization state of the circumnu- 
clear medium and/or to its covering factor. Higher quality 
data are needed to establish if only one or both of these 
parameters change in this source. By combining the infor¬ 
mation provided by the X-ray light curve and by the spectral 
analysis, we derived some constraints on the absorber loca¬ 
tion under the hypothesis of a homogeneous and a clumpy 
medium. In both cases, the absorber should be located in¬ 
side the outer edge of an extended torus and, in particular, 
under the clumpy hypothesis, it should be located close, or 
just outside, to the BLR zone. 


ACKNOWLEDGMENTS 

This work is based on observations obtained with the Swift 
satellite and it made use of data supplied by the UK Swift 
Science Data Centre at the University of Leicester. We 
thank Neil Gehrels, Boris Sbarufatti, Gianpiero Tagliaferri 
and the Swift Mission Operation Genter to make every ef¬ 
fort to get our ToO request scheduled. This research has 
made use both of the Palermo BAT Gatalogue and database 
operated at INAF - lASF Palermo and of the 70 month 
Swift-BAT Catalogue. We, moreover, thank Miguel Perez- 
Torres for the analysis of archival VLA data. Part of this 
work was supported by the European Commission Sev¬ 
enth Framework Programme (FP7/2007-2013) under grant 
agreement no. 267251 Astronomy Fellowships in Italy (As- 
troFIt). The authors acknowledge financial support from the 
Italian Ministry of Education, Universities and Research 
(PRIN2010-2011, grant no. 2010NHBSBE). Support from 
the Italian Space Agency is acknowledged (contract ASI 
INAF 1/037/12/0). Finally, we you would like to thank the 
anonymous referee for the useful and constructive comments 
which improved the quality of the paper. 

REFERENCES 

Antonucci R., 1993, ARA&A, 31, 473 
Avni Y., 1976, ApJ, 210, 642 

Ballo L. et ah, 2015, ed. by P. Caraveo, P. D’Avanzo, N. 
Cehrels and G. Tagliaferri; Proc. of “Swift: 10 years of 
discovery”, 2-5 December 2014, Rome, Italy; Proceed¬ 
ings of Science (PoS, Trieste, Italy), vol. SWIFT 10, 
PoS(SWIFT10)122 (arXiv: 1505.02593) 

Baumgartner W. H., Tueller J., Markwardt C. B., Skin¬ 
ner G. K., Barthelmy S., Mushotzky R. F., Evans P. A., 
Gehrels N., 2013, ApJS, 207, 19 
Bennert N., Jungwiert B., Komossa S., Haas M., Chini R., 
2006, A&A, 459, 55 

Bentz M. C., Peterson B. M., Netzer H., Pogge R. W., 
Vestergaard M., 2009, ApJ, 697, 160 
Bianchi S., Maiolino R., Risaliti G., 2012, Adv. Astron., 
2012, 782030 


Bohlin R. C., Savage B. D., Drake J. F., 1978, ApJ, 224, 
132 

Braito V., Reeves J. N., Gofford J., Nardini E., Porquet 

D. , Risaliti G., 2014, ApJ, 795, 87 

Burrows D. N., et ah, 2005, Space Sci. Rev., 120, 165 
Gonnolly S. D., McHardy 1. M., Dwelly T., 2014, MNRAS, 
440, 3503 

Grenshaw D. M., Kraemer S. B., George 1. M., 2003, 
ARA&A, 41, 117 

Gusumano G. et ah, 2010, A&A, 524, A64 
Elvis M. Risaliti G., Nicastro E., Miller J. M., Fiore F., 
Puccetti S., 2004, ApJ 615, L25 
Evans P. A. et ah, 2007, A&A, 469, 379 
Evans P. A. et ah, 2009, MNRAS, 397, 1177 
Goodrich R. W., 1995, ApJ, 440, 141 
Kalberla P. M. W., Burton W. B., Hartmann D., Arnal 

E. M., Bajaja E., Morras R., Poppel W. G. L., 2005, A&A, 
440, 775 

Keel W. G., 1996, ApJS, 106, 27 

Malkan M. A., Gorjian V., Tam R., 1998, ApJS, 117, 25 
Moretti, A. et al. 2005, in Siegmund O. H. W., ed., Proc. 
SPIE Conf. Ser. Vol. 5898, UV, X-ray, and Gamma- 
Ray Space Instrumentation for Astronomy XIV. SPIE, 
Bellingham, p. 360 

Munoz Marm V. M., Gonzalez Delgado R. M., Schmitt 
H. R., Cid Fernandes R., Perez E., Storchi-Bergmann T., 
Heckman T., Leitherer C., 2007, AJ, 134, 648 
Osterbrock D. E., 1989, Sky Telesc., 78, 491 
Porquet D., Reeves J. N., O’Brien P., Brinkmann W., 2004, 
A&A, 422, 85 

Reeves J., Done G., Pounds K., Terashima Y., Hayashida 
K., Anabuki N., Uchino M., Turner M., 2008, MNRAS, 
385, L108 

Risaliti G., Elvis M., Nicastro F., 2002, ApJ, 571, 234 
Risaliti G., Elvis M., Eabbiano G., Baldi A., Zezas A., 2005, 
ApJ, 623, L93 

Risaliti G., Elvis M., Fabbiano G., Baldi A., Zezas A., Sal- 
vati M., 2007, ApJ, 659, Llll 
Risaliti, G., Salvati, M., Elvis, M., Fabbiano, G., Baldi, A., 
Bianchi, S., Braito, V., Guainazzi, M., Matt, G., Miniutti, 
G., Reeves, J., Soria, R., Zezas, A., 2009, MNRAS 393, LI 
Tarter G. B., Tucker W. H., Salpeter E. E., 1969, ApJ, 156, 
943 

Trippe M. L., Crenshaw D. M., Deo R. P., Dietrich M., 
Kraemer S. B., Rafter S. E., Turner T. J., 2010, ApJ, 
725, 1749 

Whittle M., 1992, ApJS, 79, 49 


© XXXX RAS, MNRAS 000, [HS] 


